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I. INTRODUCTION 

During the last decade an increasing attention is given 
to a new class of structures called topological insulators 
(TI) with promising characteristics both in fundamen- 
tal aspects of their physics and possible applications in 
nanoelectronics, spintronics, and fabrication of new mag- 
netic, optical and information processing devices^ - — The 
principal features of TI include the presence of time- 
reversal (TR) invariance in the system where the prop- 
agating edge states may exist being localized near the 
boundary of host material and have the dispersion re- 
lation which is linear near the origin of their quasimo- 
mentum (Dirac-like structure), corresponding to the en- 
ergies belonging to the insulating gap of the bulk ma- 
terial. The spin of such states is firmly attached to 
the direction of propagation along the edge, making 
them protected against backscattering due to the TR in- 
variance which leads to effective cancelling of two scat- 
tered states with opposite possible directions of the spin 
flip which accompanies such backscattering. The exis- 
tence of such edge states have been shown in numerous 
theoretical models of TI, and also in the experiments. 
The materials included graphene^ HgTe/CdTe quan- 
tum wells^-iS bismuthjiiiH LiAuSe and KHgSb^ as 
well as general two-dimensional models of paramagnetic 
semiconductors^ the silicene; 19 ' 20 the topological modal 
semimetals.— Also a lot of studies has been devoted 
to the general properties of two- and three-dimensional 
models of TI with certain symmetries r^ -18 ' 20 " — where 
four topological invariants have been found in 3D TI in- 
stead of the single Z2 invariant in 2D TIii&i£ 

Recently a general group theoretical analysis has 
been made for the links between the geometry of the 
Bravais lattice and the properties of TI<2£ It should 
be mentioned that the symmetry arguments have al- 



ways played a significant role in classifying the sys- 
tems as trivial or topologically protected against external 
perturbationsi 6 ' 15 i 25 ~ -21 The studies of three-dimensional 
materials were mostly focused on the Bi2Se3, Bi2Te3 or 
Bi 2 Te2Seii 2 ii22i _ — where also the edge states were con- 
structed explicitly in several models of the finite-size 
geometry] 31 ' 32 It can be seen that although the features 
of the TI are very general and describe a truly novel state 
of matter, the number of different materials demonstrat- 
ing these features is currently quite limited. So, it is 
of interest to find new materials and compounds where 
possible manifestations of TI may be present, for both 
fundamental aspects and applied purposes. It is of inter- 
est also to understand which properties of the edge states 
are common for different systems, and which are special, 
and how all of them are related to the bulk quantum 
states in a specific model. 

Here we consider a model of edge states and con- 
nect their properties to the topological characteristics 
of the host material for a new candidate for the class 
of topological insulators: the two-dimensional electron 
gas (2DEG) in a material with strong spin-orbit cou- 
pling (SOC) formed at the interface of the monolayer 
of bismuth deposited on the silicon. This material is 
characterized by a giant SOC splitting which was ob- 
served experimentally also in a number of metal films 
or the combined materials of the "metal on semiconduc- 
tor" type^^r— and recently described theoretically! 38 ' 39 
Its huge spin splitting together with the hexagonal type 
of the lattice creates a certain potential of manifestation 
of the TI properties, since the spin-resolved bands may 
evolve into spin-resolved edge states, and the hexagonal 
type of the lattice is favorable for the TI to exist 

The properties of 2DEG at the Bi/Si interface have 
been studied experimentally with the help of angle- 
resolved photoemission spectroscopy (ARPES)^ ~ 38 ' 40 ~ — 



2 



applied also to other materials. It was found that this 
material represents an example of nowadays widely stud- 
ied class of materials with large (up to 0.2 .. . 0.4 cV) SOC 
spin splitting of their energy bands, which can be formed 
in various compound materials or heterostructures of the 
"metal on semiconductor" type. It is known for many 
years that that SOC plays an important role in forma- 
tion of the TI properties* 4 ^ including the localization 
effects of Rashba SOC combined with electron-electron 
interaction,— the Dirac-cone surface states in I^Sejr^ 
and B^Te^Sea-a,— PbSb 2 Te 4 or Pb 2 Bi 2 Te 2 S3^ and 
Bii-^Sbz^ the topological phase d 49 ' 50 and quantum 
spin Hall phase in honeycomb lattice^ the ultracold 
Fermi gases,— the spin Hall effect in graphene,— and 
the Kondo insulator effects . 54 ' 55 

Various materials with strong SOC have been a subject 
of intensive studies throughout recent years, including 
the structures of Bi deposited on Si-Ge superlattices,— 
the Pb on Si structure^ the trilayer Bi-Ag-Si system^ 5 - 
the structures with monolayer of Pb atoms covering 
the Ge surface,— or the Pb on Ge structures.— One 
can mention also new types of triple bulk compounds 
with strong SOC like GeBi 2 Te 4 ^ BiTel or other bis- 
muth tellurohalidesj£°-~— or recently discussed Bi^Rh^Ig 
material.— 

In the present paper we adopt the nearly-free model of 
two-dimensional bulk states in Bi/Si developed earlier— 
and applied in the extended form in our previous paper 
for the description of spin polarization, charge conduc- 
tivity and optical properties of this promising material^ 9 - 
for the calculation of ID the edge states of the electrons 
on the Bi/Si interface in the finite strip geometry. We 
obtain both the explicit form of the edge wavefunctions 
and the edge energy spectrum, calculate their spin polar- 
ization and Kubo diagonal charge conductance, and link 
the possible topological stability of their properties to the 
Zi topological invariant studied by the analysis of the 
time reversal matrix elements behavior in the Brillouin 
zonei 1 ^^ ' 15 ' 17 The results of our paper are of interest 
for expanding the knowledge of new materials with the 
topologically protected properties where the SOC plays 
a significant role, making them suitable for further ap- 
plications in spintronics as stable current-carrying and 
spin-carrying channels. 

The paper is organized as following. In Sec. II we briefly 
described the nearly free-electron (NFE) model of 2D 
bulk states at Bi/Si interface, and derive a model for 
the ID edge states for the electrons in the finite strip 
geometry. We calculate their spectrum, wavefunctions, 
and spin polarization. In Sec. Ill we apply the finitc- 
tcmpcrature and finite collision rate Kubo formula for 
calculations of the diagonal charge conductance of edge 
states as a function of the Fermi level position. In Sec. IV 
we reinforce our findings on the edge state stability by 
considering the topological band properties of the 2D 
bulk states in Bi/Si, and find the results supporting the 
presence of the TI phase. Our conclusions are given is 
Sec.V. 



II. MODEL FOR BULK STATES AND EDGE 
STATES 

A. Bulk states 

Our model for the ID edge states is based on the 2D 
NFE model for the bulk states of 2DEG formed at the 
interface of the trimer Bi/Si(lll) structure 3 ^ developed 
for the description of the spectrum near the M-point of 
the BZ, and later extended for the modeling of the elec- 
tron states in the whole BZ^S This model was compared 
with its expansion containing the anisotropic terms in the 
NFE model as well as with an empirical tight-binding 
models While the details of band structure and the 
quality of reproduction of experimental ARPES data on 
energy bands in Bi/Si vary from model to model, the sim- 
ple NFE model allows to reconstruct the main properties 
of spin split bands including the magnitude of splitting, 
the energy gap, and the spin polarization. It also has a 
major advantage of a straightforward derivation of the 
edge states in a finite strip geometry, as we shall see be- 
low. 

In this model, the Hamiltonian H = Hq + V is writ- 
ten as a sum of free-electron term Hq corresponding to 
Rashba paraboloid centered in the Gamma point of the 
hexagonal BZ, and the crystal potential energy V with 
the symmetry of the hexagonal lattice. The parameters 
of both the Hq and V are fitted as to provide the best cor- 
respondence between the model and the experimentally 
known structure of the bands near the Fermi level^ 5 - - — 
The corresponding Bloch electron wavefunctions are two- 
component spinors of the following form: 

*^*(r)=£a»*Vft? fc (r) (1) 

71 

where the Rashba eigenstates ip n u = V'k+Gn have the 
form of free electron states with the quasimomentum is 
shifted to the certain Gamma point and 

^bulk = f ±e iAr S (k y -ik x ) J ■ ( 2 ) 

The (±) sign corresponds to two eigenvalues for 
Rashba energy spectrum E(k) = h 2 k 2 /2m ± auk, and 
the material parameters are the same as in the previ- 
ous papers on Bi/Si states ; 38 ' 39 obtained for the best fit 
to the experimental ARPES data, namely, m = 0.8mo, 
a R = 1.1 cV • A and V = 0.3 cV. 

In FigUjwc plot the energy spectrum of 2DEG at the 
Bi/Si(lll) surface in the NFE model^^ Four lowest 
bands are shown with the Fermi level position reported 
experimentally 3 ^ to be between the bands 2 and 3 where 
a gap of around 0.2 eV is formed, creating the possibility 
of the edge states existence within this gap and making 
the system a possible new candidate for the topologi- 
cal insulator class. The large metallic-like values of elec- 
tron energy and SOC amplitude present for the 2DEG 
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FIG. 1: (Color online) Energy spectrum of 2DEG at the 
Bi/Si(lll) surface in the nearly free-electron model. Four 
lowest bands are shown with the Fermi level position reported 
experimentally^ to be between the bands 2 and 3 where a gap 
is formed creating the possibility of the edge states existence 
within this gap, and making the system a possible new topo- 
logical insulator. 



in this system make it promising for the consideration 
in transport and optical experiments where the disor- 
der, collision and thermal broadening prevent the ap- 
plication of conventional semiconductors. It should be 
noted that the discussed properties of the band struc- 
ture for Bi/Si 2DEG are obtained in the framework of 
one specific model with a set of parameters chosen for 
the best fit to experimental data. Thus, it may be mod- 
ified in the future when more insight will be gained on 
the properties of Bi/Si or other similar compounds. Still, 
we shall see below that the qualitative and topologically 
described features of the electron states studied within 
this model are robust against the significant variations of 
the model parameters, which is an indication of certain 
intrinsic and stable properties of the system. 



B. Edge states 

We now turn our attention to the construction of the 
model for the edge states localized at the opposite edges 
of the finite strip formed in the 2DEG. We can start 
with the strip geometry where the electrons are confined 
along the y direction in the strip —L/2 < y < L/2 and 
with conventional assumption of the hard-wall boundary 
conditions $(x,y = ±L/2) = Q^r22M^. 

First, the spectrum of edge states can be found by 
solving the eigenstate problem with the requirement of 
exponential dependence across the strip direction Oy. It 
can be done by starting from the bulk Hamiltonian and 
replacing the quasimomcntum component in the direc- 
tion of confinement by the pure imaginary variable de- 



scribing the inverse localization depth corresponding to 
the localized states exp(±Az/) which means in our case 
the substitution k y — > — iA. The eigenfunctions of this 
Hamiltonian can be constructed in the same nearly free- 
electron approximation as the bulk states (JTJ, and have 
the following form: 

* kxA (x,y) =e A ^ A (x), (3) 
Fk x A(x) = a n (k x , A)0 W fc a A(g)- (4) 

n 

The spinors 4> n k x A(x) can be obtained from ^ with 
the substitution k y —> — iA, 

0nfc x A(a;) = y= [ ±e iArg[-(fc x +nb+A)] ] ■ (5) 

The summation in ^ is over the ID lattice in the re- 
ciprocal space which corresponds not to 2D hexagonal 
but the ID simple lattice along the Ox direction with 
the period a = 2ir/b where b = 1.08 A^ 1 is the reciprocal 
space period in the k x direction. The state @ remain to 
be Bloch function along Ox with the conventional trans- 
lational property ^ kmA (x + a,y) = e lk * a <& koiA {x, y), while 
along the confinement direction the wavefunctions are 
exponents exp(±Ay). 

The wavefunction satisfying the boundary conditions 
on a single edge ^(x, y = ±L/2) = and having the spe- 
cific energy E = E(k x ,A) in the bulk gap has the form of 
linear combinations of (0]) with different exponents cor- 
responding to the given energy E = E(k x ,A): 

*kAx,y) = ^c A $ fcxA (a:,2/). (6) 

A 

If we solve the Schrddingcr equation for our model 
of 2DEG at Bi/Si interface and with the substitution 
k y —iA, the spectrum of edge states will be obtained 
as a function of two quantum numbers (k x , A). If there 
are eigenstates with energies corresponding to the gap in 
the bulk spectrum, one can be interested in them as in 
potential candidates for the edge states with topological 
protection. We have found that in our model that such 
energy levels for the edge states indeed exist in the gap 
of the bulk spectrum, and their dependence on the (fc^A) 
quantum numbers is shown in Figj^l The energy de- 
pends on two quantum numbers (k x , A), so the spectrum 
in FigHk is plotted in a 3D layout where the bulk gap 
edges are shown by two horizontal lines. The existence 
of two Dirac-like ID crossing branches of energy (pairs of 
linear dispersion curves) as functions of k x is clearly seen 
at the line k x = of the ID Brillouin zone, making these 
edge states promising for being topologically protected. 
An important feature of the edge state spectrum is the 
presence of two roots Ai^ for each value of energy for a 
given k x , i.e. the equation E(k x ,A) = E$ has two pairs 
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FIG. 2: (Color online) (a) Energy spectrum of edge states 
located in the bulk energy gap of the Bi/Si(lll) surface. The 
energy depends on two quantum numbers (k x ,A) before the 
boundary conditions are applied. The existence of two Dirac- 
like ID crossing branches of energy with the spin polarization 
S y connected to the direction of propagation is clearly seen, 
making these edge states promising for being topologically 
protected. The spin-up states are moving to the right (v x > 0) 
while the spin-down states are moving to the left (v x < 0), 
shown in the same color with corresponding branches of the 
edge spectrum, (b) Edge state energy dependence on the 
inverse localization length A for the cross-section k x — 0.05 
A . For the given position of the Fermi level in the bulk gap 
there are two roots ±Ai,2 for each edge of the strip giving two 
edge wavefunctions belonging to the corresponding branches 
of energy spectrum. 



of solutions ±Ai,2 for the left and right edge of the strip, 
respectively, as it is shown in Figl^b for the cross-section 
k x = 0.05 A . Such structure of energy eigenvalues is 
the direct consequence of the relative proximity of two 
branches of Rashba spectrum present in the basis for the 
Hamiltonian which can be also seen for the bulk spectrum 
in FigfTJ This feature allows one to construct the edge 
states satisfying the boundary conditions for the band of 
energies located in the bulk energy gap, as it is done in 
various models of edge states in TIi£ & 31 i 32 In Fig(2] we 
also plot the mean values of the only nonvanishing spin 
component S y = \ a y \ for the edge states which is 
clearly seen to be connected with their direction of mo- 
tion along the strip. The spin-up states are moving to 
the right with the group velocity v x = h~§§~ > while 



FIG. 3: (Color online) Edge states localized on the opposite 
borders of the strip (shown schematically) for the energy of 
the edge state equal to the Fermi level inside the bulk gap 
Ef — 1.5 eV and for the strip width L = 10 nm. The edge 
states are well-localized at the corresponding edge of the strip. 
For each edge there are two states propagating to the positive 
and negative directions of the Ox axis and having opposite 
spin polarizations S y shown nearby the arrow labeling the di- 
rection of movement with the color as in Fig[2] The backscat- 
tering is possible only with the spin flip, making these states 
effective for the charge transport along the edge. 

the spin-down states are moving to the left (v x < 0), 
shown in the same color with corresponding branches of 
the edge spectrum. It is known that this chiral prop- 
erty of edge states makes them topologically protected 
against backscattcring in the TR-invariant system, pro- 
viding a substantial contribution to the overall conduc- 
tance in the transport experiments^!^ 

The specific boundary conditions are applied to the 
general form of the edge state ((6]). The two wavefunc- 
tions for the given energy E = E(k x , Ai^) satisfying the 
boundary condition f = L/2) = on the left edge 

of the strip y = L/2 (when facing in the forward direction 
of the Ox axis) and decaying into the strip, has the form 
(|6]) which can be constructed explicitly by the following 
superposition of the states (jlj: 

y) = F kxAl (x) (e h ^ - eCAi-Aatf+A^ ; (7) 



^\x,y) = F kxAz (x) (e^ - e( A - A ^ +A «) (8) 

where the normalization condition is implied in 
Fk^hix). Correspondingly, the wavcfunction for the right 
edge y = —L/2 can be written as 



*<£\x,y) = F fex _ A2 (x) (e- A ^ - eCAs-AOf-A^ (10) 
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An example of the edge wavefunctions is shown in Figj3] 
for the energy of the edge state equal to the Fermi level 
inside the bulk gap E = Ep = 1.5 cV and for the strip 
width L = 10 nm. One can see that the edge states are 
well-localized at the corresponding edge of the strip on 
the length of about 1 nm. The arrows indicate the direc- 
tion of propagation and the spin S y of each state in pair 
©,© and ©,011]). The direction of the propagation of 
two chiral states on the one edge y = — L/2 in our model 
is the same as on the other edge y = L/2 which is dif- 
ferent from the most of the TI models. This difference 
can be explained due to the strong Rashba SOC in our 
system which leads, for example, to the dominant polar- 
ization S y of the states moving in the Ox direction and 
where the non-compensated total spin S y can be accumu- 
lated if the population of right-moving and left-moving 
electrons is unbalanced, for example, by external electric 
field^ which is also can be expected for the edge states 
shown in Fig(3] 

The form of spin polarization shown in Fig[2] creates 
a positive expectation about their topological stability 
for the charge transport against the scattering on non- 
magnetic impurities which do not violate the TR sym- 
metry. If we consider the backscattering, then it is clear 
from Fig[2] that the change of the propagation direction 
will lead to the spin flip, and such reversal cancels the 
reflected waves and extinguishes the backscattering?^ 
This is consistent with the arguments of the topologi- 
cal stability of such edge states as the participants of 
charge transport, which is the required property of a sys- 
tem to become a topological insulator. We shall see be- 
low in Sec. IV that our assumption about the 2DEG on 
the Bi/Si interface as a possible TI is supported further 
by the analysis of the topological properties of 2D bulk 
states. 



III. CHARGE CONDUCTANCE 

An important property of the edge states is their con- 
tribution to the total conductance of the sample which 
was considered in a number of papers on TI ^ 31 ' 32 It is 
of interest to consider its dependence on the Fermi level 
position inside the bulk gap which can be tuned exper- 
imentally. Here we perform such calculation for the ID 
edge states © taking into account the finite collision 
broadening and considering various temperatures which 
can be present in future experimental or technological 
setups with the 2DEG at the Bi/Si interface in a finite 
strip geometry. 

We start with the conventional Kubo formula^ for the 
DC conductance tensor in a system with volume f2 
obtained at finite collision rate 1/t and with the Fermi 
temperature function fi in the i-th band with energy Ei : 




FIG. 4: (Color online) The diagonal charge conductance a xx 
is plotted in units of e 2 /h as a function of the Fermi level 
position inside the bulk gap counted from the crossing point of 
two branches of the edge state dispersion curves for (a) T — 
K (solid curve), (b) T = 45 K (dashed curve), and T = 293 K 
(dotted curve) . The finite collision rate and finite temperature 
lead to the maximum values of a xx significantly below the 
ballistic value e 2 /h. When the Fermi level passes above the 
crossing point Eq of the edge state spectrum, new vacant 
states are no longer available for the electric field-induced 
transitions, and the conductance saturates. The temperature 
increase only smoothes the conductivity curves due to the 
large SOC and band gaps in the Bi/Si system being of the 
order of 0.1 eV which exceeds the scale of thermal energy. 



_ ze ^ jj - jj n it 

^ ~ n ^ Ek - E j (Ei - E,Y + (h/r) 2 ^W""*' 

(11) 

where (v v )ij are the matrix elements of the velocity 
operator v v — which contains the SOC ("anoma- 

lous") term due to the presence of Rashba SOC in the 
Hamiltonian. 

The off-diagonal (Hall) conductivity a xy vanishes since 
our system is not in a magnetic field, and is internally 
nonmagnetic. The diagonal charge conductivity a xx cal- 
culated by Eq. pTj) is plotted in FigfJ]as a function of the 
Fermi level position inside the bulk gap counted from 
the average crossing point E = Eq of two branches of 
the edge state dispersion curves (see Fig[2|). We consider 
various temperatures T = 0, 45 and 293 K, and a real- 
istic value of collision rate 1/t = 10 12 sec -1 , although 
the precise values of scattering parameters in the Bi/Si 
2DEG are yet to be determined from the experiments. 

One can see in FigQ] that the conductivity of the edge 
states even with the inclusion of the collision rate (i.e. in 
a non-ballistic regime) reaches a significant value in units 
of conductance quanta e 2 /h which makes the edge state 
in Bi/Si a promising candidate for the charge transport 
in the bulk insulating regime which has been already re- 
ported as being significant in the experiments on TI.— 
It should be noted that the maximum of the conductivity 
in Fig[J]is below the ballistic value e 2 /h. The reason is 
that we consider a non-ballistic regime with finite colli- 
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sion rate 1/r and finite temperature, both producing a 
diminishing effect on the conductivity obtained from the 
Kubo formula. Such drops of values for the edge state 
conductivity below e 2 /h have been predicted also for the 
1D2 and models of edge state where finite temper- 
atures were taken into account. 

The value of a xx rises when the Fermi level is mov- 
ing from the lower edge of the bulk gap to the crossing 
point Eq for the Dirac-like dispersion, and after pass- 
ing through this point the conductance saturates. This 
can be explained by considering the transitions from the 
lower (occupied) to the upper (vacant) states in the edge 
spectrum which all contribute to the conductance as long 
as the interband transitions dominate for the velocity 
matrix elements in (fTTj) . When the Fermi level goes 
above the crossing point in Figf2J new vacant states are 
no longer available for the dominating direct transitions 
since there is no more empty bands for the edge states 
is left above the Fermi level, so the total conductance 
tends to saturate. Since the energy spectrum shown in 
FigHk has a form of Dirac flat cone only along the k x di- 
rection and is described by a essentially different depen- 
dence along the A direction as it can be seen in Figl^b, 
the behavior of conductivity with respect to Fermi level 
variations is different from the pure Dirac spectra both 
in f D and 2D models of One can see from Fig@] 

that the temperature increase mainly smoothes of the 
conductivity curve without the major changes in its de- 
pendence due to the large SOC and large band gaps in 
the Bi/Si system of the order of 0.1 eV which exceed the 
scale of thermal energy even at room temperature, mak- 
ing the Bi/Si a promising candidate for the stable charge 
transport setups with the help of edge states in future 
room-temperature devices. 

IV. TOPOLOGICAL PROPERTIES OF BULK 
STATES 

It is known from the general theory that the stability 
of edge states is guaranteed by certain topological prop- 
erties of the bulk states. In particular, the system can 
be a TI if an integer called Zi invariant is different from 
zeroi^&iii There are several ways for calculation of this 
invariant, and here we shall use the method proposed by 
Kane and Mele& which links the Z2 index to the zeros 
of the Pfafhan for the interband matrix elements of the 
time reversal operator which has the form = ia y K for 
the spin-1/2 particles, and K is the complex conjuga- 
tion operator. If we have only two lowest bands occupied 
which is the case of our Bi/Si 2DEG (sec FigJTJ) , then the 
Pfaffian Pi^k) is equal to the single off-diagonal matrix 
element between the Bloch functions 111,2 (k) in the occu- 
pied bands 1 and 2, 

Pi, 2 (k) = (ui(k) 1 e I «a(k)>. (12) 
The topological considerations provide a convenient 




FIG. 5: (Color online) Absolute value of the Pfaffian |T2]) 
in the hexagonal BZ for the 2DEG on the Bi/Si interface 
for different values of bulk band parameters, (a) Vo = 0.3 
eV, an = 1.1 eV-A; (b) Vo = 0.6 eV, a R = 1.1 eV ■ A; 
(c) Vo = 0.3 eV, a R = 0.6 eV ■ A. The Pfaffian has three 
pairs of zeros at the corners of the BZ where the visible zeros 
are shown as black circles while at TR-invariant T and M 
points shown as blue and green circles the value | Pi t 2 \= 1. 
These properties indicate that the Z2 invariant is odd, and 
the topological insulator phase is present. 
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form of using the definition (| 12[) for the challenging task 
of finding new material demonstrating the properties of 
TI. Namely, if for the hexagonal BZ the k-depcndent 
function -P(k) has pairs of zeros in the corners of the 
BZ (or, depending on the overall symmetry, on the lines 
inside the BZ), then the system demonstrates the prop- 
erties of TIi^ i 11 ! 15 There is an extensive discussion of the 
Zi invariant properties related to the TI including an- 
other definition of it where the matrix elements of TR 
operator are calculated between the states with oppo- 
site wavevectors k and — k and the TI is determined by 
its properties not in the entire BZ but at discrete set of 
"time reversal invariant points" like T or M points. The 
detailed discussion and all relevant mathematical connec- 
tions between different approaches on the calculation of 
Z2 invariant can be found in the literature* ^ 15 ' 17 ^ 7 ' 66 " — 

In Fig(5]a, we plot the absolute value of the Pfaffian (|T2|) 
in the hexagonal BZ for the 2DEG on the Bi/Si interface 
for the same basic set of model parameters as were used 
for calculations of the bulk spectrum in FigfTJ In order to 
see the possible changes in the Z2 index with the varia- 
tion of the system parameters, in FigJ3]3,c we plot | P\ : 2 \ 
for two other sets of parameters: in FigjSjD the ampli- 
tude of the periodic potential is increased compared to 
the initial case shown in FigJT] Vo = 0.6 eV and the am- 
plitude of Rashba SOC is the same, = 1.1 eV • A. In 
FigJSJi the periodic potential amplitude is the same as in 
FigEk, Vq = 0.3 eV, but the Rashba coupling amplitude 
is decreased, cur = 0.6 eV • A. It is clearly seen for all 
cases that the Pfaffian has zeros in the corners of the BZ 
where the visible zeros are shown as black circles which 
border is shown schematically, while | Pi. 2 |= 1 in the 
TR-invariant T and M points shown as blue and green 
circles, respectively. There are three pairs of zeros for 
I -Pi, 2 I which indicates that the Z2 invariant is odd, thus 
classifying the 2DEG at the Bi/Si interface as a topolog- 
ical insulator with the protected edge states ,L2i£iiii. 

One can see in FigJSJ that the variations of material 
parameters do not change significantly the topological 
properties of the Pfaffians which all have the same quali- 
tative features with | Pi 2 |= 1 at the TR-invariant T and 
M points and with three pairs of zeros for | Pi 2 | in the 
corners of the BZ. The depth of the parameter variation 
present in three parts of FigJS] is rather big and reaches 



50% which covers a wide range of possible experimental 
realization of the 2DEG at the Bi/Si interface. Still, the 
absolute values of Pfaffians shown in FigJS] look very sim- 
ilar to each other which indicates their qualitative topo- 
logical nature being the key for discovering new examples 
of topological insulators. The method of mutual analysis 
of chiral edge states and topological bulk properties used 
in our calculations can be applied to other materials and 
structures. 

V. CONCLUSIONS 

We have derived a model for the one-dimensional edge 
states for the electrons on the bismuth on silicon inter- 
face in the finite strip geometry. Based on the bulk nearly 
free-electron model, their energy dispersion was obtained 
inside the bulk gap being linear in the quasimomentum. 
The spin polarization of edge states is linked to the direc- 
tion of propagation along the given edge which provides 
topological stability of these chiral modes. The diagonal 
charge conductance of the edge states was calculated as 
a function of the Fermi energy for finite collision broad- 
ening and finite temperatures, and has shown a specific 
behavior with respect to the Fermi level position rela- 
tive to the edge state dispersion curves which may help 
in their experimental detection. The topological stabil- 
ity of edge states was confirmed by the structure of the 
interband matrix element for the time reversal operator 
which was shown to be stable against the large variations 
of material parameters. The results of the paper may 
be of interest both for the development of the topologi- 
cal insulator theory by providing a novel example of the 
material belonging to this class, and also for the develop- 
ment of new spintronics and nanoelectronics devices with 
stable transport and operating at room temperature. 
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